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Quantum stochastic operator equations are derived for inter-branch exciton and polariton pro- 
cesses caused by acoustic phonon scattering. The use of a fully quantum model combined with these 
recently developed techniques predicts the presence of "stimulated scattering" terms, and provides a 
sound basis for understanding the basis of the approximations used in generating the equations. The 
theory is applied to a model motivated by recent experiments where a stronger photoluminescence 
signal from the high energy microcavity polariton polariton is observed for low excitation powers, 
and also to an experimental scheme designed to show evidence of the stimulated scattering effect. 



I. INTRODUCTION 



The development of stochastic operator and wavefunc- 
tion approaches to dissipative processes caused much in- 
terest in the field of quantum optics. My aim here is to 
present a derivation of equations for inter-branch phonon 
processes using one of these techniques, that of quan- 
tum stochastic differential equations (QSDE's) jij. The 
derivation of these equations relies on similar approxi- 
mations to the well known master equation techniques 
jj], 0], and so has the advantage of being a systematic 
derivation in which the approximations are clearly stated 
and whose effects are generally well understood. Master 
equation techniques from a "quantum optical" viewpoint 
have already been applied to both excitonic optical bista- 
bility [|| and the semiconductor laser B . The stochastic 
wavefunction approach based on an "unravelling" of the 
master equation |q has been applied to both electron- 
phonon scattering [3| and exciton-phonon scattering [Q. 
This exciton-phonon work was for a model involving a 
single quantum-dot exciton level, and the emphasis was 
on non-Markovian effects. 

Here the QSDE technique is used to treat the phonon 
interaction as a coupling between two exciton or polari- 
ton branches and an infinite reservoir of phonon modes. 
Polaritons are formed in a semiconductor quantum mi- 
crocavity, structures in which quantum wells are grown 
inside an optical cavity. The light field is confined by two 
Bragg mirrors, and the quantum wells create confined ex- 
citon states. With appropriate design, the photons in the 
optical cavity and the excitons couple together creating 
a quasi particle called a cavity polariton. The QSDE's 
for the polariton system include factors accounting for its 
partly photonic nature. In addition, they contain extra 
"stimulated scattering" terms, also predicted by Pau et 
al jllj] , whose presence would not have been predicted by 
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a phcmoncnonological approach to modelling the system. 

This paper was motivated by photoluminescence (PL) 
measurements of microcavity polaritons at low tempera- 
tures (1.7 - 50K) and a range of input laser intensity (0.4 
- 10 Wcm" 2 ) |h]. Most investigations of these cavity po- 
laritons use reflectivity measurements, with only a few 
considering the other optical emission properties. Con- 
trary to expectations, it was the higher energy polariton 
that dominated the PL at the lowest laser power densi- 
ties. The experimental results are due to the thermali- 
sation of the free carriers down into the upper polariton 
branch and the dominance of the cavity decay rate over 
the scattering between the polariton branches. Free car- 
rier scattering becomes insignificant at low carrier densi- 
ties and LO phonon processes do not occur because their 
energies are greater than the polariton splitting. As a re- 
sult, the photon intensity (PL signal) emitted from decay 
of the upper polariton is strong, while the signal from 
the lower polariton is restricted by the small influx of 
excitation due to the slow inter-branch acoustic phonon 
scattering. The equations for the theoretical model are 
solved in an appropriate limit to predict the severity of 
this bottleneck. The theory is also extended to a sys- 
tem where the upper polariton is coherently excited by a 
laser tuned to the polariton energy, and predicts a sim- 
ilar bottleneck effect. The results depend on the size of 
the photonic component of the polariton as well as its 
excitonic part. 

Finally, the stimulated scattering term is investigated 
by constructing a model experimental system that could 
be used to measure its effects. This model also predicts 
the strength of the effect, which is proportional to the 
exciton- phonon coupling, the exciton- photon mixing, 
and the intensity of the stimulating field. 



Note added (March 2002) 

Stimulated emission of bulk polaritons was experimen- 
tally observed in II- VI sample in the late 1970's [18|, 
and microcavity polaritons in QW samples in 1997 [19|. 
Stimulated emission terms (i.e. a factor of N + 1) were 
inserted by hand into the theory of Tassone et. al. |2(|] 
on the bottleneck effects in microcavity polaritons, and 
then with more theoretical justification in by Tassone and 
Yamamoto [^lj which included discussion of stimulated 
scattering. In these two papers the interest is restricted 
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to generating a rate equation model, and the papers pay 
little theoretical attention to all the steps in a rigorous 
derivation such as I have attempted here, and certainly 
quantum noise terms are not considered at all. 

The bottleneck effects and stimulated scattering have 
become of much recent interest - see for example some 
of the references brought to my attention by J. J. Baum- 
berg |^2|, |23|, Q . It turns out the model I used when 
doing these early calculations at University of Sheffield 
in 1996 is far too simplistic, and so is not relevant to 
any likely experimental work. However, the theoretical 
method used to generate both the QSDE's and model is 
still valid, and it would be a straightforward matter to 
recast the model to a more appropriate one and do some 
useful calculations. 



II. THE PHONON INTERACTION 

In this section I introduce a model Hamiltonian for 
the phonon scattering process. This does not specify 
the type of phonon explicitly, they could be either op- 
tical or acoustic (either longitudinal or transverse), as 
the phonons are treated as a reservoir of boson modes 
with particular properties, properties that are necessary 
for the QSDE technique to be valid. As long as the type 
of phonons considered have these properties, the results 
given here are applicable. 

It is also relevant to consider here how we describe the 
exciton and cavity modes. For low excitation strengths 
the exciton can be described as a boson. This is be- 
cause the resulting low density of excitons means that 
any phase space filling || or exciton-exciton scattering 
processes are unlikely to occur. In addition, there is a 
continuum of exciton modes for different in-plane mo- 
menta. I will refer to this as an 'exciton branch', and the 
energies of these are given by the dispersion relations for 
the exciton. Similarly, a microcavity photon mode can 
also have a range of k values, and can also be called a 
'photon branch'. 

Consider an interaction Hamiltonian for quantum well 
excitons interacting with phonons. The two dimensional 
nature of each exciton branch is allowed for by a set of 
exciton modes indexed by the in-plane wave vector k, 
an these modes are associated with the operators e,(A;), 
ej(fe). The phonon modes are indexed by in-plane wave 
vector k as well as their growth direction wave vector 
q, with the operators b(k,q), b'(k,q). The coupling be- 
tween exciton and phonons is denoted x- I n the inter- 
action picture, the operators all have an oscillatory time 
dependence given by their mode energies, and so hcua^d 
type terms do not appear in the Hamiltonian. This time 
dependence is not explicitly included to avoid cluttering 
the notation. The general form of the interaction Hamil- 
tonian for first order scattering involving any two exciton 
branches and one phonon branch is 



ij k k' k" q 

[*(*)+ e+(fc)] [e](k') + e,(k') 
~b\k",q) + b(k",q) 



l (k)+e\(-k) 



(2.1) 



but 



Often the terms are written like 

this is simply a re-ordering of the same summation. The 
i, j indexes span the exciton branches, the k summations 
are over the in-plane wave vector of the excitons, and the 
q sum is over the wave vector of the phonon in the growth 
direction. Note that with only minor modifications to 
the summations and indices the same expression can be 
used to describe phonon scattering in bulk or even 1 or 
dimensional structures. 

The first simplification we can make is to apply con- 
servation of in-plane crystal momentum k, which is not 
conserved because of the lack of translational symme- 
try in the growth direction. This reduces the number of 
summations in the above formula from three to two, and 
results in 



Ha = EEEE^( fc ' fc - Afc '9) 

ij k Ak q 



[ei(fc)et(fc- Ak)P(Ak,q) 

+e\(k)ej(k- Ak)b(Ak,q) 
+e ] l {~k)e](k - Ak)tf(Ak,q) 

+ii(-k)ej(k - Ak)b(Ak, q) 



(2.2) 



It is normal in quantum optics to make a rotating wave 
approximation to remove non-resonant (non energy con- 
serving) terms because these off resonant terms oscillate 
at about twice the optical frequency. This is very much 
faster than any likely dynamical process in the system, 
and so their net effect is assumed to average to zero, in 
a kind of "coarse graining" approximation. In this sys- 
tem involving phonons, it is not obvious that this type of 
approximation can be made. Both the excitons and the 
phonons have a continuum of energies, and because the 
phonons can have frequencies near zero, the oscillations 
in the off resonant terms can easily be relatively slow, so 
the rotating wave approximation (RWA) cannot always 
be applied. A situation like this can be dealt with using 
non-rotating wave techniques (eg. juj) 

In the case of optical phonon scattering a RWA is 
possible. This is because optical phonons have a min- 
imum frequency (about 10 13 Hz), and the off resonant 
terms would be oscillating at about twice that rate. Also, 
acoustic phonon scattering could be treated between two 
sufficiently separated exciton branches, as again the off 
resonant terms would oscillate at twice the frequency of 
the separation of the two excitons. Alternatively, if for 
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some reason the exciton to phonon coupling vanished suf- 
ficiently fast for small phonon energies the exciton could 
only couple to relatively energetic phonons for which a 
RWA might be reasonable. Assuming well separated 
branches or optical phonons, the oscillations in the off 
resonant terms will occur on much shorter timescales 
than any other dynamical process in the system. This 
means we can apply the RWA, so the Hamiltonian be- 
comes 



Ha = j2J2J2J2^ k ~ Ak ^~) 

i^j k Ak q 

[ei(k)e](k - Akp(Ak,q) 
+e](k)e j (k- Ak)b(Ak,q) 

A. The phonon QSDE 



(2.3) 



Following Gardiner jTj, we can write down the Ito 
quantum stochastic differential equations (QSDE's) for 
the above model. In this approach we treat the phonon 
modes as a reservoir or 'heat bath' with which the exci- 
tons can interact. At first sight, this QSDE method does 
not seem directly applicable to this problem. This is be- 
cause there is only one reservoir of phonons in the system, 
whereas the standard application of the theory would as- 
sume that each scattering between excitons of different 
wave number is coupled to an independent reservoir of 
phonons. For this reason some correlations might per- 
sist (or perhaps form) between different exciton modes 
because they interact with this common reservoir. We 
might reasonably assume this to be an insignificant ef- 
fect because the phonon reservoir is much bigger than 
the range of exciton states - the phonon reservoir is three 
dimensional, whereas the exciton states only fill two di- 
mensions. 

For example, compare this with the standard quantum 
optics case, where a single system mode interacts with a 
one dimensional continuum of reservoir modes - there 
are an infinite number of reservoir modes for the sys- 
tem mode to decay into. Similarly, the exciton phonon 
scattering model presented here also has sufficiently large 
reservoirs in which to "lose" any quantum correlations. 
Although plausible enough, this argument is not suffi- 
cient. Leaving the detailed calculations to a later work 
here I just make the extra assumption that any such cor- 
relations will decay away, and so the use of QSDE's here 
is valid. 

There are three properties this reservoir must have for 
the QSDE approach to be valid. Firstly, there must be 
a smooth, dense spectrum of phonons with which each 
exciton mode can interact. For acoustic phonons, which 
have an energy proportional to the magnitude of their 
wave vector, such a spectrum exists. For LO phonons, 
the phonon modes only exist above a certain energy. This 
can cause problems due to the effect of the edge of the 



LO phonon band, but if the inter-exciton transition is 
sufficiently larger than the lowest LO phonon energy, the 
scattering will not see the phonon band edge. 

Secondly, the coupling between system and reservoir 
must be linear in the reservoir operators, in agreement 
with the first order interaction Hamiltonian used above. 

Thirdly, the coupling constants between the reservoir 
and system must be smooth functions of frequency, which 
is holds for deformation potential coupling to phonons 
|Q. This leads to the additional and related approxi- 
mation usually called the First Markov Approximation. 
This requires the coupling parameter in the Hamiltonian 
between the excitons and the reservoir of phonons to be 
independent of frequency in the range of interest. This 
"range of interest" is where energy is very nearly con- 
served - where the phonon carries away (or supplies) the 
energy lost (or gained) by the conversion of the exciton. 
Note that the coupling is assumed to remain roughly con- 
stant over the width of the resonance. This is also a kind 
of "coarse graining" , as a flat frequency response implies 
a strongly peaked time response - so the coupling appears 
"flat" if the width of the time response is much shorter 
than the important dynamical processes occuring in the 
system. Also, the state of the reservoir must be unaf- 
fected by the interaction. In this case, the energy lost 
by the excitons in the form of phonons must not change 
the state of the reservoir - for example, by increasing its 
temperature. This is guaranteed for low excitation levels 
by the relative sizes of the system and reservoir. 

The procedure has allowed us to eliminate the phonon 
modes from the model. It has given us a system in which 
the phonon scattering causes a loss of energy from the 
excitons and added a delta correlated "quantum noise" . 
This noise is the effect of the quantum uncertainty of the 
phonon reservoir feeding into the exciton system, and is 
described by new quantum noise operators dBp, dBp. 

The Ito QSDE for some arbitrary system operator S 
affected by phonon scattering between two different ex- 
citon branches is 



ds = ^2Y1 Xij ^ k,k ~ Ak ^ 

iytj Ak 2 

{[TVp(Afc) + I] [2eJ(fe)ey(fe - Ak)Sei(k)e](k - Ak) 
-Se\{k)e 3 {k - Ak)ei{k)e]{k - Ak) 
-e\{k)e 3 {k - Afc)e i (jb)et(jfc - Ak)S 
+N P (Ak) \2e t (k)e](k - Ak)Se\{k)e 3 {k - Ak) 
-Sei(k)e](k - Ak)4{k)ej{k - Ak) 
-ei{k)e]{k - Ak)e\{k)e 3 {k - Ak)S } 

st 



dt 



E E k - Ak) S, e]{k)e 3 (k - Ak) 

i.j Ak 



dB T 
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E E yXij (k, k - Ak)dB P S, e t (k)e](k - Ak) 



i.j Ak 



(2.4) 



This operator S could be any operator of interest - per- 
haps one of the mode operators, the intensity operator 
or even some kind of correlation function. The opera- 
tors dBp, dBp are quantum white noise operators, and 
Np(Ak) is the average number of thermal phonons in 
each phonon mode. The summation over q has dropped 
out because the resonance condition allows us to calcu- 
late it in terms of the k terms and the dispersion curves. 
The noise operators dBp, dB P are denned by their cor- 
relations, which are 



dBp(t)dB P {t'\ 
dBp(t)dBp(t'] 
dB P (t)dB P (t'] 



dBp(t)UB P {t')* 
= (N P + l)5(t-t')dt, 
= N P 5(t-t')dt. 



= 0, 



(2.5) 



These operators cause "quantum white noise" , and cor- 
respond to a reservoir in which the number of thermal 
quanta is constant per unit bandwidth. This is not the 
same as a quantum thermal noise, which would have 
a Boltzmann distribution for its thermal quanta. Al- 
though using quantum white noise is not exact, it is 
a very convenient and widely used approximation. In 
particular, it would not be a good approximation for 
the low wavevector part of an acoustic phonon reservoir, 
where the large relative change in energy between nearby 
wavevectors means a proper thermal ensemble should be 
used Q. However, since we have already introduced a ro- 
tating wave approximation bases on a large energy sepa- 
ration between exciton branches, and hence are not in the 
low energy part of the phonon spectrum, using quantum 
white noise does not cause any extra restriction. 



B. Intensity and amplitude operator equations 

To get a clearer picture of what the phonon interaction 
does, I will give some examples relevant to the calcula- 
tions in this paper. In the quasi-PL and coherent excita- 
tion models presented in section VI, equations for excita- 
tion number and mode amplitude are used; more specif- 
ically, these are derived from the QSDE's for the mode 
intensity and mode amplitude operators. For simplicity, 
I write £i(k) as e+, Bj(k — Ak) as e_, and Xij{k,k — Ak) 
as X- The intensity for the (+) mode is just /+ = e^_e+, 
that for the (— ) mode is 7_ = ete_. The modes de- 
scribed by e+ and e_ are different. 

The QSDE equations for the intensities of the (+) and 
(— ) modes due to the phonon scattering are 



dl+ = x -I+U 



1 - N P 1+ - I. 



y/\ e\_e-dB P + dB P e+el 



dl- = X 



N P I+- J_ 



dt 



+t/x e\e-dB P + dB P e + el 



(2.6) 



The noise terms cannot be rewritten in a form depend- 
ing solely on I±, so I leave them in terms of e+, e„. 
However, we could define new operators J\ = e^e_ and 
J2 = e + eL, and from them we get extra equations which 
allow us to remove the e type terms. Note how the terms 
in the pair of equations balance each other to preserve 
excitation number, a consequence of the fact that the 
scattering transfers excitation between exciton branches 
while emitting (or absorbing) phonons. This contrasts 
with the usual loss models that use reservoirs which de- 
stroy excitations by removing them from the system. 

The QSDE equations for the amplitudes of the (+) and 
(— ) modes due to the phonon scattering are 



de+ 



1 - N P e+ 



dt + yfxz-dBp, 



ele+e_ — Np£- 



dt + y/xdB P e + . (2.7) 



dt 



These equations have an interesting feature - the pres- 
ence of a "stimulated scattering" or "incoherent stimu- 
lated emission" term. These show up as the I+I- and 

Np (l_ — in the intensity equations. A naive con- 
sideration of the phonon scattering process would not 
suggest such terms, as the process is incoherent and is 
caused by loss. Hence it is natural to expect only simple 
loss terms or population transfer terms. The new terms 
are predicted by this quantum model because the two 
exciton modes are coupled together by the phonon inter- 
action, and so their bosonic nature causes these stimu- 
lated effects. However, because the phonon scattering is 
Markovian, the stimulated emission loses its phase mem- 
ory and becomes "stimulated scattering" . The appear- 
ance of these terms shows clearly the advantages of work- 
ing from a more rigorous derivation, as opposed to a phe- 
nomenological approach based on simply writing down 
rate equations. They were also predicted by Pau et al 
JTT|] who started from the Heisenberg equations for this 
system. 



III. AN IDEAL POLARITON 

A polariton consists of an exciton coupled to a cav- 
ity mode. A full model would include a description of 
the different types of exciton, their excited states, as well 
as all the cavity modes and the couplings between these. 
However, in practise we know that excitons will only cou- 
ple to a cavity mode if they have matching wavevectors. 
Also, for low excitation strengths only the lowest energy 
exciton is significant. This means I can consider an "ideal 
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polariton" , consisting of an idealised quantum well exci- 
ton, an idealised cavity mode, and the exciton to pho- 
ton coupling term between the two. Further, I use the 
Coulomb gauge, the dipole approximation, and the ro- 
tating wave approximation to express the coupling term 
in its simplest form. 

The result couples an exciton to a cavity mode, where 
both are modelled by quantised harmonic oscillators. I 
write et(fc), e(k) as the operators for the exciton mode, 
and a^(fc), a(k) as the operators for the field in the cavity 
mode. The interaction Hamiltonian is 



e(k) = c + (k)p + (k) +c_(fc)j5_(fc) 

a(k) = d + (k)p+(k) + d-(k)p^(k). (3.5) 

These last equations are the most useful for the work in 
this paper. They allow us to take an interaction involving 
only the exciton or photon modes and rewrite it in terms 
of the eigenstates of the coupled system, ie. in terms of 
the polariton operators. 



Hp = hA(k) [e\k)a{k) + e(k)tf(k)] (3.1) 

The parameter A(k) is the coupling between the cavity 
mode and the exciton. It is in units of frequency, and is 
equal to half the splitting between the pair of polaritons 
at resonance. I will denote the exciton energy as Kl{k), 
and the cavity mode energy as hw(k). This two mode 
system of the coupled exciton and cavity mode can be 
diagonalised, and then be written in terms of new po- 
lariton operators, which describe the two new polariton 
modes that can be used in place of the original exciton 
and photon modes. The two eigenvalues (or frequencies) 
of the eigenmodes (the polaritons) are 



2e±(fc) = fi(fc) +u(k) 

±^[n(k)-uj(k)} 2 + 4A(k) 2 . (3.2) 

The polaritons are characterised by a label ± denot- 
ing the "upper" or "lower" (in energy) polariton branch, 
and their wavenumber k. If I introduce the operators 
p±(k), p±(k) to describe the polaritons, the eigenvectors 
of the system describe how the cavity and exciton op- 
erators combine to form the polariton operators. The 
polariton operators for the coupled exciton - photon sys- 
tem are pj_(fc), p±(k). The exciton and cavity modes 
have coupled together to form two new polariton modes. 

The polariton operators are related to the exciton and 
cavity operators by the following expressions 



p+(*0 



c+(fc)e(fc) 
c_(A)e(fc) 



■d+(k)a(k) 
■d-{k)a(k) 



(3.3) 



where 



c+(fc) = -d-(k) = — =yj I + sm [a(k)] 
v2 



-1 



d + (k) = c_(fc) = —Vl-sinMfc)] 



tan [a(k)] 



[O(fc) - u(k)] 
2A(k) 



(3.4) 



Note that the equations for the polariton operators in 
terms of the exciton and cavity operators can be inverted. 
The results are 



IV. COHERENT DRIVING, LOSSES, AND 
THERMAL NOISE 



Of course it is not only the phonon scattering which 
is main topic of this paper that needs to be considered. 
Any realistic model needs to include losses due to exciton 
recombination and the imperfect cavity mirrors. Also, 
sources of excitons or polaritons need to be included - 
these can originate either from the combination of hot 
free carriers to generate a thermal population of exci- 
tons, or by coherent excitation from a laser beam. To 
do this I use the quantum stochastic differential equa- 
tion (QSDE) for a polariton with coherent driving, linear 
losses and thermal noise. When describing a realistic po- 
lariton, a sum over in-plane wave number k is needed to 
account for the dispersion. However, since each k mode 
is independent in these interactions, I do not need to in- 
clude the sum over k states here. The standard driving 
and loss terms for a boson mode described by operators 
/, P have an interaction Hamiltonian like 



H. 



D 



(4.1) 



The first term describes a coherent driving process, 
which can be imagined as being due to a classical driving 
field, but is actually a coherent state driving term. The 
second term couples the mode to a reservoir consisting of 
an infinite number of boson modes ft, f j. The loss from 
the system is modelled by irreversible transfer of excita- 
tion from the system into this reservoir. The amount of 
loss is described by a decay rate 7 that is related to the 
7j', but includes the sum over the density of states of the 
reservoir. 

In a polariton system the excitons are not the eigen- 
states. The eigenstates are the polaritons, and so I need 
to replace the exciton operators with their e quiv alent in 
terms of polariton operators, defined by eqn fl3.5[) . I will 
now use f± to represent c± or d± , depending on whether I 
am considering interaction with the cavity mode (/ = d) 
or the exciton (/ = e). In polariton operators, the inter- 
action Hamiltonian can be rewritten 



Hn = 



{[e*f + ]p + (t)-[e*f + rpl(t)} 
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+i{[e*f-]p-(t)-[e*f-}*pUt)} 

i 

+j2[f-p-wt+f-pUm 



(4.2) 



The QSDE for an arbitrary system operator A affected 
by this Hamiltonian is just pj. 



dt 

, 7 



{e*f + }p + -{e*f + }*p%,A 



dt 



[N + 1] [/;/+] \2p\A P+ - Ap%p+ - pIp+a 



+ | [N] [/;/+] [2p+Apl - Ap+pl-p+p^A 



V7 



A flP 



t 



dB + yfydBi A, f+p 



(4.3) 



The [...] represent the p_ terms, which are similar to 
the p+ ones, and N is the thermal occupation number 
of the reservoir modes. Different sets of reservoir modes 
will have different thermal populations according to a 
Boltzmann weight factor. A reservoir coupled to a mode 
of frequency u will have a thermal occupation iV pro- 
portional to exp [— huj/kBT]. At optical frequencies and 
room temperature, this factor is so small that there are 
usually no significant thermal effects. The quantum noise 
operators dB, dW arc defined in the same way as those 
in eqn (2.7). 



A. Intensity and amplitude equations 

To get a clear picture of what these interactions do, 
I will calculate some specific examples of these QSDE's. 
In the quasi-PL and coherent excitation models presented 
in section VI, equations for excitation number and mode 
amplitude are used; more specifically, these are derived 
from the QSDE's for the mode intensity and mode ampli- 
tude operators. First I consider the equation of motion 
for polaritons with the exciton (or photon) affected by 
losses and thermal noise, but with no coherent driving 
term. The QSDE for how the intensity I± = p±p± of 
the mode changes due to the losses and the associated 
thermal noise is 



of the equation, with the quantum noise terms averaged 
to zero. Replacing the expectation value with a number 

I± = (l±\ the equation can be rewritten as 



N] 



(4.5) 



Now I consider the effect of a coherent driving term as 
well as losses. The QSDE for how the operator p± of the 
mode changes due to the losses and the associated noise 
is 



dp± 



n 



f±f±-jP± 



dt + y/yfedB. (4.6) 



A complex number equation for the mode amplitude 
can be produced from this QSDE easily. We can take the 
expectation values of the equation, and assume the mode 
is in a coherent state \a). The quantum noise term then 
averages to zero, and we get an equation in terms of the 
coherent amplitude a 

While this approach might seem fine if the mode is 
in a coherent state, what if it is not? Fortunately, the 
basic idea can still be used. Any quantum state can be 
expanded over a basis of coherent states Q , b y using ei- 
ther the Glauber-Sudarshan P representation |IJ, |f5), or 
by the subsequent generalisations |l6). Here the system- 
atic procedure is to convert the QSDE into the equivalent 
master equation form and then apply the representations 
in the standard way. In practise this means that the state 
of the mode is not specified by a single complex number 
a, but by an appropriately distributed ensemble, each 
member of which follows the equation 



-a = f ± e-f ±f± -a. 



(4.7) 



In these coherent state representations the non ther- 
mal part of the quantum noise terms vanish, as they 
are absorbed into the coherent state basis. The effect 
of thermal noise on the system is to spread the ensemble 
out - this does not affect its average amplitude, but will 
change its intensity by some amount. If thermal noise 
can be neglected, and the system starts or is being co- 
herently driven, then it can be adequately described by a 
single complex number. However, if there are additional 
interactions, such as strong phonon scattering, then the 
system might well depart strongly from a coherent state. 



dt 



I± = 7/±/± 



-I±+N 



dt 



v/7 \fip±dB + f±d&p ± 



(4.4) 



Rate equations can be produced from this QSDE easily. 
All we need to do is to neglect the quantum noise terms, 
and replace the intensity operators with numbers. This 
is done more rigorously by taking the expectation values 



V. THE POLARITON-PHONON INTERACTION 

In section II I described an interaction Hamiltonian for 
excitons scattering off phonons. To convert it into one 
relevant to a polariton system, I employ the same proce- 
dure as for the coherent driving and decay interactions. 
Since the eigenstates of the system are the polaritons, 
and so I need to replace the exciton operators with their 
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equi valen t in terms of polariton operators, as defined by 
eqn (3.5). A single exciton and cavity mode branch gives 
us a two polariton system, so the interaction Hamiltonian 
in the polariton operators is 

k Ak q 

{[d+(k)p + (k;t)+d-(k)p-(k;t)] 
d* + {k- Ak)p%(k- Afe;t) 

+d*_(k- Ak)pl(k- Ak;t)] P(Ak;q) 



Acoustic phonon; scattering 



+ 



d\{k)p\{k;t) 



dt(k)pl{k;t) 
[d + (k - Ak)p + (k - Ak;t) 
+d-(k - Ak)p- (k - Ak; t)} b(Ak; q)\ . (5.1) 

When this is expanded out, there are four types of 
important term. Two are intra- polariton terms, involv- 
ing phonon scattering along the upper or lower polariton 
branches, and the other two are inter- polariton scatter- 
ings. In addition to these, there are non-resonant terms 
which do not conserve energy such as those that create 
or destroy two polaritons. 

How do all these processes fit together? First consider 
the dispersion curves for the polaritons and the acoustic 
phonons. This is a 2D problem - the polariton is con- 
fined in the growth direction by the microcavity DBR 
mirrors, and the quantum well. Because both polaritons 
only have in-plane freedom, only phonons that conserve 
in-plane momentum can be emitted. The dispersion re- 
lations for the polaritons are bowl-like surfaces, and the 
phonon dispersion is like a solid cone. Acoustic phonons 
carry very little energy for their momentum, so this cone 
is very shallow compared to the other dispersions. The 
'cone' is solid because the phonons can have a non zero 
growth direction momentum - the internal parts of the 
cone correspond to when some of the phonon momentum 
is in the growth direction. 

Figure [l] is a schematic diagram showing the disper- 
sion curves and allowed acoustic phonon scattering pro- 
cesses. The intra-branch process mainly transfers popu- 
lation from large k to small k along the same branch, 
although the opposite can also occur in the presence 
of thermal phonons. The inter-branch processes trans- 
fers population between the two branches, predominantly 
from the upper (+) to the lower (— ) branch. Again, 
the reverse process is possible in the presence of ther- 
mal phonons. However, note that the upwards thermal 
processes use a phonon dispersion with its apex on the 
lower initial state. Similarly, figure ^| shows the dispersion 
curves and allowed optical phonon scattering processes. 
This follows the same principle as the acoustic case, but 
with the optical phonon dispersion curve and so includes 
a minimum phonon energy. 

Here I am primarily concerned with the downward 
inter-branch process, from the upper to the lower branch 




FIG. 1: The acoustic phonon scattering process. The thin 
curves represent the dispersion relations defining the exciton 
or polariton branches. The triangular shaded area gives the 
allowed energies and in-plane momenta of the phonons. The 
area is shaded because the phonons can have a growth di- 
rection momenta which has been projected out. The overlap 
of the branches and the shaded area is the set of possible fi- 
nal exciton or polariton states given that the initial state was 
at the dot at the apex of the phonon triangle. (A) Intra- 
polariton phonon scattering. The emission of a phonon can 
transfer the polariton from its initial state to another part of 
the polariton branch. (B) Inter-polariton phonon scattering. 
The emission of a phonon can transfer the polariton from its 
initial state to part of the lower polariton branch. 



with the production of a phonon. This means a RWA 
can be justified for sufficiently large branch separation as 
the non resonant terms will oscillate at a rate given by 
twice this - much faster than the speed of the important 
system dynamics. The resonant upper to lower phonon 
branch interaction term is then contained in 



6 =EEE^ fc ^- Afc ^) 

k Ak q 

\d+{k)p+(k; t)d*_(k - Ak)pl(k - Ak; 0(Ak; q) 

+d* + {k)pl(k;t)d_(k - Afc)p_(fc - Ak;t)b(Ak;q) 

(5.2) 

Note that this is the same as the exciton phonon 
Hamiltonian if we replace d+(k)p + (k) with e(k) and 
d*_(k- Ak)pl(k- Ak) with e(k - Ak). 



A. The polariton-phonon QSDE 

The justification for the validity of the QSDE method 
for this system is identical to that used in the initial 
exciton scattering section, so I do not repeat it here. The 
only thing that has changed is the introduction of the 
coefficients for the component of exciton in the polariton, 
and the particular allowed scattering transitions between 
branches. The Ito QSDE for a general system operator S 
in a two polariton system undergoing (+) to (— ) phonon 
scattering is 
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LO phonon scattering 



noise oper ator s dBp, dB P are defined in the same way 
as in eqn (2.7). 




k=0 



FIG. 2: The LO phonon scattering process. The thin curves 
represent the dispersion relations defining the exciton or po- 
lariton branches. The roughly parabolic shaded area gives the 
allowed energies and in-plane momenta of the LO phonons. 
The area is shaded because the phonons can have a growth 
direction momenta which has been projected out. The top of 
the parabola is displaced downwards from the dot denoting 
the initial position of the exciton or polariton by the mini- 
mum energy of an LO phonon. The overlap of the branches 
and the shaded area is the set of possible final exciton or po- 
lariton states given that the initial state was above the apex 
of the phonon parabola. (A) and (B) show the intra and 
inter-polariton scattering. 



Ak 



X {k,k-Ak) 



d+(k)d*_(k - Ak)\ 



{(N P (Ak) + 1) 



dt 



2p\{k)p-{k - Ak)Sp + (k)pL(k - Ak) 
-p%(k)p-(k - Ak)p+(k)pl(k - Ak)S 
-Sp\(k)p-(k - Ak)p + (k)pl(k - Ak) 

-Np(Ak) \2p + (k)pl(k - Ak)Sp%{k)p-{k - Ak) 
-p+(k)pl(k - Ak)pl(k)p-(k - Ak)S 
-Sp+(k)pl(k - Ak)pl(k)p-(k - Ak) dt} 
- Vx(k,k-Ak) [d* + {k)d- (k - Ak)] 



Ah 



\s,p\(k)p-(k- Ak) 
+ Vx(k,k~Ak) [d+(k)d- (k - Akf 



Ah 



d& \s,p+(k)pUk- Afc) 



(5.3) 



The operators dBp, dB p are quantum noise operators 
which replace the infinity of phonon operators b and , 
and Np(Ak) is the average number of thermal phonons in 
each phonon mode. The summation over q has dropped 
out because the resonance condition allows us to calculate 
it in terms of the fc terms and the dispersion curves. The 



B. Intensity and amplitude operator equations 

Similarly to subsection II. B, I calculate some relevant 
examples of these QSDE's, and convert them into rate 
and complex number equations. For simplicity, I write 
p+ (k) as p+ , and f>- (k — Ak) as p- . The intensity for the 
(+) mode is just I + = p\p+, that for the (— ) mode is 
/_ = pLp- ■ The polariton coefficients I will also simplify, 
writing d+(k) as d+, and d_(/c — Afc) as d_, and x(k, k — 
Ak) as simply \- 

The QSDE equations for the intensity of the (+) and 
(— ) modes due to the phonon scattering are 



dl+ = x\d+d*_ 



dl- 



X \d+d 



iJL 



1 + N P /_ - L 



dt 



^/X d + d_p\_p-dBp + d* + d-dB P p + p 

2 



1 



N P /_ 



dt 



+Vx d%d-p\_p-dBp + d + d*_dBpp + pL\ . (5.4) 



As in the exciton case, the noise terms cannot be 
rewritten in a form depending solely on I±, so I leave 
them in terms of the p+, Again the terms in the 
pair of equations balance each other to preserve excita- 
tion number. Using the same procedure as in IV. A, these 
have a rate equation form 



It 1 ' 



X \d+d_r (/_ + 1) + n p (/_ - /+)] 



= x\d + d*_\ 2 {I + (I_ + l)-N P (!_-/+)] .(5.5) 

Similarly, the QSDE equations for the amplitude of the 
(+) and (— ) modes due to the phonon scattering are 



dp^ 



X\d+d*\ -p + 



pip- + 1 + N F 



dt 



'Vxd* + d-p-dB Pl 



dp- = X \d-\-d*_ I p- 



p\p+ - N P 



dt 



-y/xd+d*_dB P p^ 



(5.6) 



Using the procedure referred to in section IV. A, these 
have a complex number form for the coherent amplitude, 
which is 



da j 



X \d+d 



* i 2 



(ula- + lj - Npaq 



dt 



2 X d* + d-Npala-dW A 



8 



Dr.Paul.Kinslcr@physics.org 

EPRINT Quantum stochastic theory of phonon scatt... http://www.kinsler.org/physics/ 



da- 



/Xd*+d. 


-(1 + 


2Np)a+a~ 


i 2 




-l_Q;_ — Npoc 


2 x d* + d 


-(i + 


Np)a f + a+(l 


lxd*+d„ 




2Np)a + a- 



dt 



r. (5.7) 

These polariton phonon scattering equations have the 
same features as the exciton equations derived in section 
II. B. This means the same comments apply here too - no- 
tably those about stimulated scattering. The new terms 
dW are white noise increments, dW a and dW E are un- 
corrected real white noises, and dWi is a complex white 
noise, where (i = {A, B,2}) 



(dW l (t)dW j {t')*) = 5 l3 5{t - t'). 



VI. POLARITON PHONON MODELS 



(5.8) 



Here I consider two models one of a 

photoluminescence-like process and the other in- 
volving coherent excitation. Both models involve a pair 
of polariton branches that result from the coupling of a 
single exciton and a cavity. The polaritons are described 
by boson operators p±(k), p±(k), indexed by in-plane 
wave vector. The splitting between the higher (+) and 
lower (— ) energy polaritons is assumed to be less than 
the minimum optical phonon energy, so only acoustic 
phonons need to be included. The coupling to acoustic 
phonons is much weaker than that to optical phonons, 
and this causes the bottleneck effect. 

The photon component of the polaritons can decay due 
to cavity losses, which I model by coupling to optical 
reservoirs. These represent the field modes outside the 
cavity. The cavity decay rate is 7i(fc), with thermal pho- 
ton number N^k), and quantum noise operators dB^k), 
dB L (k). The exciton component of polaritons on either 
branch can recombine and emit photons, which I model 
by coupling each exciton mode to independent photon 
reservoirs. This recombination rate is 75 (fc), with ther- 
mal photon number N E (k), and quantum noise operators 
dB E (k), dB^k). 

The acoustic phonon process couples the (+) modes to 
all allowed (— ) modes. The presence of thermal phonons 
is described using a thermal phonon number Np(k). Ex- 
periments done by Fainstein et al fl7f| for Raman scatter- 
ing but in an otherwise similar situation show no power 
dependence in the spectra. This indicates that the effects 
of thermal phonons are likely to be negligible, so we can 
set Np(k) = 0. Only the downwards (+) — > (— ) inter- 
polariton phonon scattering process is included here. 

I will simplify both models with some physical con- 
siderations. Both the recombination reservoirs and the 
cavity decay reservoirs have very low thermal excitations 
(Ni,(k),NE{k) w 0), because the exciton and cavity 



mode energies are much greater than the thermal en- 
ergy. Further, the recombination decay r ) E {k) is typi- 
cally a much slower process than the cavity decay rate 
lE{k) *C 7z,(fc), and so can usually be ignored when the 
polariton has appreciable fractions of both cavity mode 
and exciton. 

In the following equations, I will suppress the k ar- 
gument when k = for brevity, so that p + = p + (0), 
c± = c±(0), d± = d±(0), 7l = 7l(0), 1e = 7b(0), 
N L = N L (0), N E = N E (0), and I also write x(0, fc) as 
X(k). 



A. Photoluminescence 

The photoluminescence model is one in which the po- 
lariton is pumped as the free carriers excited by the driv- 
ing field can combine in concert with the cavity mode to 
form polaritons. Exciton formation in PL is a compli- 
cated process, with free carriers combining to form hot 
excitons (with k ^> 0), and excited excitons (eg. those 
with extra orbital momentum). These decay by LO or 
acoustic phonon emission or other scattering processes 
into the "cold" exciton states that can decay radiatively, 
and can therefore be seen optically. The situation with 
polaritons is similar. 

A full description of this process is beyond the scope 
of this paper, so here I use an approximate model. I 
suppose that the recombination processes preferentially 
fill the upper polariton branch, something which could 
occur because the upper branch has higher energy than 
the lower branch for the same k. This agrees wth the 
results that originally motivated this paper Js) in which 
the upper branch is preferentially populated. Further, 
I assume the statistical properties of the newly created 
polaritons are influenced most strongly by the thermal 
nature of the free carriers they originated from. We sup- 
pose that the predominant effect of the optical pump- 
ing is that each mode (k value) of the upper polariton 
branch is driven by an independent thermal population, 
modelled by coupling to a set of reservoir modes which 
cause a weak loss 7c(fc) ~ 0, but have a high thermal 
occupancy Nc{k) so there is significant thermal driving 
strength Mc(k) = ^fc{k)Nc{k). The free carrier decay 
rate is small because the polaritons are unlikely to be 
scattered into free carriers. 

I make the further simplifying assumption that the fi- 
nite k higher polariton states decay rapidly down into 
the k = state due to optical and acoustic phonon scat- 
tering. As a result only the k = upper polariton state 
is relevant, and I couple a single k — upper polariton 
to a reservoir described by the parameters 7c, Mc, and 
the quantum noise operators dBp, dB p . 

The complete equations for I + = l+(0), I-(k) can 
then be constructed. Firstly, however, I can simplify 
them with the physical considerations given above. For 
the (+) mode the sources of thermal noise are swamped 
by the free carrier component, since Mc = 7c-^c 3> 
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photoluminescence 




I_(k) 



k=0 



FIG. 3: The quasi- photoluminesence model. The thin 
curves are the polariton branches. The triangular shaded area 
gives the allowed energies and in-plane momenta of the acous- 
tic phonons. The overlap of these two is the set of possible 
final states given that the initial state was at k = 0. The thick 
arrows represent excitation cascading down the upper branch 
by primarily LO phonon scattering to give a thermal popu- 
lation near the k = upper state. The thick lower curve is 
the allowed final states after emission of an acoustic phonon. 
The dashed arrows are the light emitted by decay from the 
polaritons through the cavity mode. 



"iL{k)NL{k)^E{k)NE{k). Also, the free carrier decay 
rate is much smaller than the cavity decay rate 7c <C 
7l (k) . These simplifications hold true for all cases where 
the polariton has appreciable fractions of both cavity 
mode and exciton. 

The phonon scattering itself is very weak, so that 
7L ^ X)fc x(k) \d+dL(k)\ . In the typical situation where 
J_(/c) remains small because the lower polariton is only 
populated by this scattering process, this means that the 
effect of the phonon scattering terms in the equation for 
I + is negligible. 

Figure || shows a schematic diagram of the model. The 
QSDE's for it are 



dU 



dl-(k) 



{- lL \c + \ 2 i++M c }dt 

+^lEc + p\dB L (Q) 

+^Ec + dB L (0)^p +1 

{- lL (k)\c_(k)\ 2 I + 



(6.1) 



+x(k)\d + d*_(k)\ z i+ 



L(k) + l}} 



dt 



+ VlL(k)c*_(k)p f _{k)dB L {k) 

+ ^ lL {k)c„{k)dB L (k) ] p+(k) 

+ ^Mtyd* + d-(k)plp-(k)dB P (0) 

+ Vx(k)d + d*_{k)dBl(0)p + pl(k). (6.2) 

Since in this application we do not need the detailed 
quantum mechanical noise correlations, we can approxi- 
mate them by taking the expectation values and factoris- 
ing the different intensity moments (see section IV. A). 



Replacing the expectation value with a number I± 
I± \ , I can describe the system using rate equations 



d I 



-7i(0) \c + \ I+ + M c 
- lL {k)\c„(k)\ 2 1_{k) 
+x(k)\d + d_(k)\ 2 I + {I_(k) + l}} -(6.3) 
The exciton modes have steady state populations of 



J_(fc) 



Mr 



\C+\ 1L 

X(k) \d+d- 



\c-{k)\lL{k) 



(6.4) 



where the equation for /_ (k) has been further simpli- 
fied by reusing the fact that the phonon scattering is a 
weak process, with |c_(fc)| 2 7i(fc) 3> J2k x(k) \d+d-(k)\ 2 
and /+ small. Each polariton mode emits photons at a 
rate given by the product of its decay rate |c±(fc)| 2 7x,(fc) 
and the occupation I±. For the upper (+) mode this is 
just Mc- Using the (+) mode population, the (— ) modes 
have populations of 



J_(fc) = 



X (k)\d+d.(k)\ 2 M c 
|c + | 2 7 L|c-(fc)| 2 7L(fc)' 



(6.5) 



The ratio of cavity decay intensities between all (— ) 
modes and the k = (+) mode is 



R = 



E 



X (k)\d + d-(k)f 

\c+\ 2 lL 



(6.6) 



From this expression we can see that the relative 
strengths of the upper and lower polariton emission peaks 
is controlled largely by the strength of the phonon scat- 
tering between the two branches. As a result, the upper 
peak is larger than the lower peak. In contrast, a model 
that coupled the free carriers to both branches would 
populate them according to their Boltzmann factors, and 
would have the lower exciton peak as the larger because 
of its greater thermal population. 



B. Coherent excitation 

Now I will now consider a system where the higher en- 
ergy polariton (+) is coherently excited by a laser beam 
tuned to k = 0. The laser field is assumed to be strong 
enough so as to be unaffected by the energy lost in cre- 
ating polaritons. This does not necessarily imply a 'high 
power' experiment, since even low intensity beams have 
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coherent excitation 




I_(k) 



k=0 



FIG. 4: The coherent excitation model. The thin curves 
are the polariton branches. The triangular shaded area gives 
the allowed energies and in-plane momenta of the acoustic 
phonons. The overlap of these two is the set of possible final 
states given that the initial state was at k = 0. The left 
hand side dashed arrow represents the coherent excitation 
E — |e+| 2 of the upper k — state, and the two right hand 
side arows are the light emitted by decay from the polaritons 
through the cavity mode. 



many photons, and the coupling between field and polari- 
tons is weak. The rate that new polaritons are created is 

€+. 

The (+) mode is dominated by its coherent driving, 
and so I write down the equation for the evolution of its 
amplitude. In contrast, the (— ) mode is populated solely 
through the incoherent phonon scattering, and so has no 
coherent amplitude. As a result I use an equation for its 
intensity, which is non zero. 

The complete equations for p + = p+(0), 1 ~(k) can 
then be constructed. Firstly, however, I can sim- 
plify them with the physical considerations given above 
by removing the weakest decay terms and the small 
thermal contributions. The phonon scattering itself is 
very weak. Because of this, we also have 7l(0) ^> 

Sfe x(k) ■ In the typical situation where /_ (fc) 

remains small because the lower polariton is only popu- 
lated by this scattering process, this means that the effect 
of the phonon scattering terms in the equation for p + is 
negligible. 

Figure [| shows a diagram of the model. The QSDE's 
for it are 



dp+ 
dL(k) 



dt + ^Ic* + dB L {Q) 



{-j L (k)cUk)c-(k)L 
+ X (k)\d+d*_(k)\ 2 plp + [L(k) + l]} 
+VlL(k)c*_(k)pUk)dB L {k) 



dt 



WlL(k)c- (k)dB L (k)^p+ (k) 

+^Mk)d* + d-(k)p{p^(k)dB P (0) 

+ ^X~(kjd + d*_(k)dBU0)p+pHk). (6.7) 



Since we do not need the quantum mechanical noise 
correlations, we can approximate them by taking the 
expectation values and factorising the different inten- 
sity moments (see section IV. A). Using J_ = (l—^ an d 
q + = I can rewrite the equations as complex am- 

plitude and rate equations 



!'-<*> 



1L | ,2 



+ ■ 



- 7i (fc)Mfc)| 2 /_(fc) 

X(k)\d + d^(k)\ 2 \a + \ 2 ll-(k) + l] .(6.8) 



The exciton modes have steady state values of 



a+ 
I-(k) 



2c, 



il ic+r 



X(fc)|d + rf-(fc)| 2 |a+| 2 
\c-(k)\ 2 iL{k) 



(6.9) 



where the equation for /_ (k) has been further simpli- 
fied by reusing the fact that the phonon scattering is a 
weak process, with |c„(fc)| 2 7i(fc) » J2k x(fy \d+d-(k)\ 2 
and small \a\ . Each polariton mode emits photons at a 
rate given by the product of its decay rate \ c±\ 7i(fc) and 
or i_ . For the upper (+) mode this 
Using the (+) mode population, 
) modes have populations of 



the occupation |a+ 
is just |c + (fc)| 2 7i |a+| 2 



the (- 



J_(fc) = 



X (fc)|d+cL(fc)| 2 |e + | 2 
|c + | 2 7L|c-(fc)| 2 7L(fc)' 



(6.10) 



The ratio of cavity decay intensities between all (— ) 
modes and the k = (+) mode is 



R 



E 



X (k)\d+d-(k)\' 
|c+| 2 7i 



(6.11) 



From this expression we can see that the relative 
strengths of the upper and lower exciton emission peaks 
is controlled largely by the strength of the phonon scat- 
tering between the two exciton branches. As a result, 
the upper exciton peak is larger than the lower exciton 
peak. Note that in the limits considered here, this result 
is identical to that for the PL system. 



VII. STIMULATED SCATTERING 
EXPERIMENT 

The stimulated scattering is perhaps the most inter- 
esting effect to be predicted by this quantum model of 
phonon scattering. Given this, how would it be possible 
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stimulated scattering 




k=0 



FIG. 5: The stimulated scattering experiment. The thin 
curves are the polariton branches. The triangular shaded area 
gives the allowed energies and in-plane momenta of the acous- 
tic phonons. The overlap of these two is the set of possible 
final states given that the initial state was at k = 0. The left 
hand side dashed arrow represents the coherent excitation of 
the polariton states. The upper arrow is a continuous driv- 
ing field E\ = |e+| 2 exciting the upper k = state, and the 
lower is the pulse train that periodically excites the chosen 
lower state at k = ko- The two right hand side arows are the 
light emitted by decay from the polaritons through the cavity 
mode. The upper level emission has a series of dips caused by 
the stimulated scattering induced by the periodic excitation 
of the lower polariton level. 



to measure its effects experimentally? Since the size of 
the effect is small due to the weak exciton-phonon cou- 
pling, it will be difficult to separate the small depletion 
due to stimulated scattering (the signal) from the optical 
decay from the polaritons (the background). 

Consider a simple polariton system with two branches. 
There are cavity losses, coherent driving fields tuned 
to the polariton energies, and inter-branch coupling by 
acoustic phonon processes. Recombination is neglected 
as the cavity decay is generally much larger. The upper 
branch is driven by a CW coherent field at its k = 
energy. This means that in the absence of any thermal 
phonons, there are no phonon processes along the upper 
branch to complicate matters. For simplicity I also only 
include a single k value on the lower branch k = kg, in 
an approximation which is justified at the end of this 
section. The lower level is driven periodically by very 
short laser pulses that are well separated by a time in- 
terval t, with r 3> I/71, 1/72. Between each pulse, the 
system returns back to its "relaxed state" , which is that 
for an undriven lower level. Note that this relaxed state 
is the same as that worked out in the coherently pumped 
bottleneck calculation in section VLB. Figure |] shows a 
diagram of the model. 

Choosing a kg 7^ angle on the lower branch means 
we can separate the output signal from the upper and 
lower branches by detection angle. This is important 
because the phonon effect is very small, and could be 
easily masked by the addition of a strong background due 
to the driving fields. Also, by varying the lower branch ko 
that we select, we can probe how the the phonon coupling 
strength varies with transverse wavevector k. 



The coherent driving terms mean that is is most con- 
venient to use the c- number amplitude equations. The 
equations are (Np = 0) 



dt 



7+ X 1 * -1 \ 

e+ - ya + - —a+ (a_a_ + 1) 



dt 



+iy/x'a+a-W 2 , 



+ ^/x'a* + a + W B + i^x'a+a-W; (7.1) 



The equations can be explicitly converted into those 
for a polariton system by setting a+ — * a+(0), 
a_ -> a_(fc ), 7+ = 7lc+(0)*c + (0), 7_ = 
7 L c_(fc )*c_(fco), e + = e + c +(0)*, e_ = e_c_(fc )*, x' = 
X(0, — k o )d + (0)d-(ko)* . It is also possible to do this kind 
of experiment in a two branch exciton system by using 
angled and appropriately tuned laser beams for direct ex- 
citation of the desired point on the exciton branch. Al- 
though in a PL setup a steady source of excitation might 
be supplied by the combination of free carriers into po- 
laritons, this is neither k specific, nor can it be rapidly 
modulated, nor can the individual branches be selectively 
excited. 

I now expand equations ( |7.l| ) about their relaxed state, 
and write 



a+{t) = a +0 + /3 + (t), 
ou(t) = a_ + ft(*), 



(7.2) 



where a±o are the 'relaxed states' of the two levels. 
The depletion ft of the upper level can be treated as 
a small perturbation, since it will be small compared to 
its level of excitation a+o- The relaxed state is like that 
for the bottleneck calculations so it has a_o = because 
there is no coherent driving and the excitation due to 
phonon scattering is incoherent. The intensity 7_o = 
a*_ a-o ~ x' l a +o| 2 /t+ would be non zero because of the 
W2 noise terms. Because the excitation of the lower level 
is very small so the change in its excitation ft c ann ot be 
treated as a perturba tion . Inserting equations (7.2) into 
the equations above (7.1) gives us 



1 



dt' 



7+ [-ft 



r_ 

7+ 
x' 



ft (a*_ a- + fta_ + cft ft + ft ft) 



a+o (ft a-o + "-oft + ft ft) 

7+ 



at 2 7+ 



, X 7+ / * 
1 (a +0 a +0 



7+ 7- 

+fta +0 + a+oft + ftft)] P 



(7.3) 



To do the calculation rigorously I rescale the time by 
7 + so that everything could be expanded in powers of 
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the small parameter g = x' H+- I have not defined new 
scaled variables, but just collected the factor together in 
an appropriate way. 

The truncation to first order in x' H+ is easy if 
7-/7+ ~ 1, and 0*_ 0- ~ 1. The condition on /3_ still 
allows the lower level driving pulses to be non pertur- 
bative, but ensures that the depletion they cause in the 
upper level is perturbative. Further, since the depletion 
depends on the coupling, I assume j3 + is at most of or- 
der x' M+- This assumption is self consistent with the 
resulting solutions. 

Note that these conditions are needed for a simple the- 
oretical solution for this problem, but the character of the 
results for more general conditions is unchanged. Trun- 
cating the equations to first order in x' ll+ involves ne- 
glecting any terms including two of x'/7+! 0+i a -o- The 
equations then can be easily simplified down to 



7i 
" 2 



0+ - — a+oPl 0- 
7+ 

i- 



(7.4) 



I have introduced an effective decay rate j'_ — 7_ — 
2x'a* +Q a +0 . 

Now I have set up the model, I can investigate the effect 
of the lower level driving pulses. Imagine one of the e_ 
pulses arrives at t = 0. For simplicity, I assume that the 
width of the pulse is much shorter that 1 /7_ , so it can be 
treated as a delta function. Laplace transform methods 
then quickly give the simple exponential solution for 0- 
as 




FIG. 6: Graphs of (a) a — h and (b) a_ as functions of time 
for the model without the quantum noise terms. The spectra 

for a h (c) show the sidebands caused by the stimulated 

scattering. 




FIG. 7: Graphs of (a) a h and (b) a_ as functions of 

time for the model with the quantum noise terms. The spectra 

for a h (c) show the sidebands caused by the stimulated 

scattering. 



0- (t) = a- exp (-7^/2) . (7.5) 
The equation for + is now completely defined, and is 



dt 



0+ = -?±0 + - ^ a+0 \a_\ 2 eM-"f'-t) (7.6) 



This can also be solved by Laplace transform methods, 
to give the pair of solutions 



7 

/?+ 
7 



x'a+o |a-T 
27 

27- -7+^0 



exp (-7+4/2) [1 - exp(-Tt)] , 



Xa+o \a-\' 



t exp (-7+4/2) : 



= 2i_ - 7+ = 0. 



(7.7) 



So the result is a dip in the excitation of the upper 
level, with a duration of about 1 /7+ and a relative depth 
proportional to x' \ a -\ ~ note that this is consistent with 
the assumptions that + was of order x' ll+ which was 
used in the truncation to first order. The size of the 



dips is proportional to the pulse intensity (and not its 
amplitude). A train of these input pulses exciting the 
lower level will stimulate phonon scattering in bursts and 
produce a corresponding series of dips in the upper level 
amplitude. These dips will show up in the optical decay 
from the upper mode, and could be detected by a peak 
in the output spectrum at a frequency given by their 
repetition rate. A single dip would be hard to detect, 
because it would be easily masked by noise. These dips 
show up in simulations of the equations - figures |^ and |7] 
show the behaviour of the outputs as a function of time 
for the both the simple noiseless case as well as the case 
with quantum noise included. 

Lastly, it is necessary to consider the effect of the other 
phonon scattering processes in a real system, since up to 
now we have restricted it to a simple two mode model, 
when in reality there will be scattering to other parts of 
the lower branch as well, and so the excitation in the 
chosen lower level k value will depend on re-scattered 
excitation as well as the driving pulses and the scatter- 
ing from the upper level. However, these processes will 
merely provide a nearly constant background - mostly 
due to the excitation cascading down along the lower 
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branch through k — ko down to k = 0. The only way 
they will change the signal because of the pulse train is 
through second order processes. The effect of the reduced 
excitation has to scatter down into the lower level (one 
factor of x')> then back up (another factor of x')- These 
types of processes already occur in the simple model, but 
were neglected due to their small size, and so will have 
negligible effect on the measurements. 

Light could also scatter from the pulse train directly 
into the k = detector. The strength of this could be 
measured by turning off the CW upper level driving field 
and retaining the pulse train. Also, any instability in the 
driving laser at a frequency about that of the pulse rep- 
etition should be reduced as much as possible to avoid 
contaminating the depletion signal. This should be min- 
imised by stabilising the laser, or by adjusting the repe- 
tition rate to fall in a quiet part or the laser spectrum. 

VIII. CONCLUSION 

This paper has presented a derivation of equations 
for phonon processes using the technique of quantum 
stochastic differential equations (QSDE's) Jj|. The pro- 
cedure had the advantage of being a systematic deriva- 
tion in which the approximations are clearly stated and 
whose effects are generally well understood. The tech- 
nique was used to treat the phonon interaction as a cou- 
pling between the two polariton branches and an infinite 
reservoir of phonon modes, and properly accounted for 



the partly photonic nature of the polariton. The basic 
theoretical models of experimental situations used here 
in conjunction with the weak nature of acoustic phonon 
scattering predicts a "bottleneck" effect like that seen in 
recent low power polariton experiments [|| . In addition, 
the rigourous quantum approach to the phonon scatter- 
ing interaction revealed some new "stimulated scatter- 
ing" terms whose existance would not have been obvious 
from a more phenomenonological approach. 

Further experimental investigations of this effect might 
take one of two paths. Firstly, the detuning of exci- 
ton and photon might be varied to check the change in 
strength of the bottleneck as the fractions of exciton and 
photon in the polaritons vary. The change in strength 
should mean that the ratio of the intensity peak areas 
should change in a corresponding way. The second path, 
and the one with more potential is to explore the use of 
coherent excitation. With appropriate tuning of the driv- 
ing field and perhaps even control of its angle to the per- 
pendicular, specific parts of the polariton branch could 
be excited and the resulting angular dependence of the 
luminescence measured. 

An experiment was proposed to try to observe the 
"stimulated scattering" predicted by the quantum model 
used in this paper. However, whether it is possible to 
separate this effect out from the "background" of other 
processes is debatable. Theoretical avenues would involve 
treating the intra-branch phonon scattering properly, and 
doing more detailed modelling of particular experiments. 
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